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SUMMARY 


The  objective  of  this  program  is  to  investigate  new 
a nal y t i c a  1/ nume r i c a 1  techniques  for  microstrip  antennas  and  arrays  of 
noncanonical  patch  geometries.  Antenna  performance  improvements,  such  as 
frequency  band  broadening  and  dual-frequency  operations,  are  the  ultimate 
goal  in  this  investigation. 

In  the  course  of  this  research,  a  computer  program  for  the  radiation 
of  an  arbitrary  Hertzian  dipole  in  a  general  stratified  multilayer  medium 
was  developed.  Although  the  original  intention  was  for  application  to 
microstrip  antennas,  this  algorithm  was  later  found  to  be  applicable  to  a 
broad  class  of  integrated  circuits,  geophysical  problems,  etc. 

The  analysis  of  microstrip  antennas  was  focused  on  the  development  of 
a  model  of  equivalent  magnetic  current  along  the  periphery  of  the 
microstrip  patch.  This  is  a  new  technique  derived  from  a  similar  approach 
in  dealing  with  microstrip  circuits  of  arbitrary  shapes.  Although 
numerical  convergence  problems  presently  exist,  this  approach  appears  very 
attractive  for  noncanonical  microstrip  antennas.  Continued  research  in 
this  approach  is  highly  recommended. 

Research  in  microstrip  arrays  was  conducted  for  diode-switched 
microstrip  panels  sandwiched  by  dielectric  layers.  A  moment  method 
approach  was  taken  with  the  aperture  expanded  in  terms  of  exponential  basis 
functions  and  a  single  mode  in  the  diode-gap  region.  The  computer  program 
developed  in  this  research  was  checked  against  known  data  for  the 
reflection  and  transmission  of  sandwiched  screen  structures. 

All  of  the  three  techniques  investigated  are  promising  approaches  to 
treat  microstrip  antennas  and  arrays.  It  is  recommended  that  substantial 
efforts  be  devoted  to  extend  and  approve  the  techniques.  These  future 
efforts  should  lead  to  progress  in  microstrip  antenn3  technology  as  well  as 
geophysics  and  integrated  circuits. 
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SECTION  I 


INTRODUCTION 


The  requirement  for  Low-profile  antennas  with  low  fabrication  costs 
has  led  to  the  recent  interest  in  microstrip  antennas.  However,  microstrip 
antennas  have  been  found  to  be  limited  in  their  performance,  especially 
regarding  their  typically  narrow  bandwidth  of  one  to  six  percent.  A 
theoretical  investigation  of  noncanonical  antennas  appears  promising  in 
developing  antennas  of  higher  performances,  as  suggested  in  several 
experimental  studies.  In  this  study,  techniques  for  analyzing  microstrip 
antennas  and  arrays  were  investigated  in  order  to  develop  computer 
algorithms  for  noncanonical  patch  shapes. 

Three  major  accomplishments  are  reported  in  the  following  section!..  A 
general  numerical  method  for  analyzing  the  radiation  of  an  arbitrary 
Hertzian  dipole  in  multi-layer  stratified  media  was  developed.  Roth  the 
number  of  layers  and  their  thicknesses  and  complex  dielectric  constants  can 
be  arbitrarily  designated.  The  Hertzian  dipole  can  be  vertical  or 
horizontal,  electric  or  magnetic.  The  computational  method  was  developed 
as  a  basis  for  the  microstrip  antenna  analysis.  However,  its  application 
was  later  found  to  be  quite  broad,  including  geophysical  problems  and 
integrated  circuits. 

The  method  of  analysis  for  the  microstrip  antenna  developed  in  this 
research  is  based  on  the  representation  of  fields  in  terms  of  an  equivalent 
horizontal  magnetic  current  along  the  periphery  of  the  microstrip  patch. 
This  method  is  particularly  attractive  for  patches  of  arbitrary  shapes.  An 
integral  equation  in  terms  of  equivalent  peripheral  magnetic  current  was 
derived  and  solved  by  the  method  of  moments.  Although  the  numerical  data 
do  not  show  the  expected  convergence  behavior,  the  basic  formulation  and 
computer  algorithm  represent  a  highly  promising  approach  which  should 
provide  a  powerful  tool  for  the  analysis  of  noncanonical  microstrip 
antennas  when  fully  developed. 

Analysis  of  microstrip  diode  switched  arrays  is  a  difficult  prob'  si 
for  which  little  research  has  been  performed.  In  this  research,  a  moment 
method  computi  r  program  was  developed  for  a  simple  d i ode- swi t c he. I 
monolithic  slot  (or  strip)  array.  This  program  has  been  tested  again;! 
data  for  a  screen  structure  which  have  been  presented  in  the  literature. 


These  three  major  accomplishments  represent  progress  in  geophysical 
science  stratified  media  problems  and  integrated  circuits  as  well  as  tie- 
intended  microstrip  antenna  and  array  analyses.  All  of  the  analyses 
involve  new  and  useful  numerical  techniques  which  need  further  improvement 
for  full  f rui t i on . 
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APPLICATION  OF  HERTZIAN  DIPOI-E  ANALYSIS  TO 
MICROSTRIP  ANTENNAS 


A.  Introduction 

A  general  mierostrip  antenna  is  depicted  ir,  figure  I.  The  feeds, 
shown  in  dotted  lines,  can  be  either  a  probe  feed  or  an  end  feed.  In 
general,  the  end  feed  is  preferred  because  it  is  simple  to  construct.  Both 
the  geometry  of  the  microstrip  patch  and  the  location  and  type  of  feed  are 
parameters  determining  the  characteristics  of  she  mierostrip  antenna. 
Current  mierostrip  antenna  technology,  which  has  been  extensively  reviewed 
in  an  IEEE  Transactions  on  Antennas  and  Propagation  Special  Is.su..:  [  1  j  and 
two  reference  books  [  2,3],  does  not  have  analytical  methods  suitable  for 
patches  with  noncanonical  shapes. 

The  study  of  mierostrip  antenna  (excluding  arrays)  technology  has 
progressed  from  rectangular  patches  [4,5]  to  circular  [  6],  elliptical 
[7,8],  and  triangular  [5  ]  patches,  etc.  With  a  new  geometry,  there  is 
usually  a  new  set  of  performance  characteristics  that  represent  an 
improvement  over  existing  designs  and  may  even  find  different  new 
applications.  However,  as  summarized  by  Carver  and  Mink  [  9  I,  the 
mierostrip  antenna  has  typical  bandwidths  from  one  to  six  percent.  Larger 
bandwidths,  which  may  be  achieved  by  increased  substrate  thickness  or  patch 
size,  are  often  impractical.  Kernweis  and  Mcllverma  110,11)  discovered 
that  mierostrip  antennas  could  operate  at  dual  frequencies  if  radial 
conducting  strips  were  added  to  the  patch. 

The  method  of  analysis  for  mierostrip  antennas  is  typically  a 
technique  highly  dependent  on  the  geometry  of  the  patch.  In  fact,  it  deals 
as  a  rule  mierostrip  patches  with  canonical  shapes  conforming  to  the 
several  existing  coordinate  systems.  Some  efforts  were  made  with  limited 
success  to  apply  t  he  wire-grid  method  j  17]  and  the  finite-element  tech...  que 
[13]  for  put. 7:  ■;  <>:  generalized  shapes. 

In  the  present  analysis,  we  begin  with  the  analysis  of  the  fields 
to  an  orbit  i.iry  Hertzian  dipole  in  an  arhitiat'y  stratified  medium.  This 
approach  is  in  essence  similar  to  the  Green's  functions  formulation  such  as 
that  of  Cavalc ante ,  et  al  l 1 4 ! .  The  numerical  analysis  techniques  involved 
in  this  approach  hr.  been  reported  in  the  Interim  Report  1  15  |  nnd.ue  being 
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formally  published  ll6  J. 

Two  methods  were  studied  in  this  research  to  analyze  microstrip 
antennas,  using  ttie  Hertzian  dipole  radiation  as  a  building  block.  One  is 
to  replace  the  conducting  microstrip  patch  by  a  sheet  of  equivalent  surface 
current.  The  other  is  to  formulate  the  problem  in  terms  of  the  peripheral 
electrical  fields  of  the  patch.  These  are  discussed  as  follows. 

B.  Method  of  Replacing  the  Microstrip  Patch  by  an  Equivalent  Current 

Sheet 

The  conducting  microstrip  patch  can  be  replaced  by  a  sheet  of  current 
density  £5,  as  shown  in  Figure  9.  This  is  accomplished  by  using  the 
surface  equivalence  theorem  on  the  patch  surfaces  and  then  allowing  the 
thickness  of  the  patch  to  approach  zero.  The  vertical  J  symbolizes  a  probe 
feed.  For  an  end-fed  patch  the  excitation  current  should  be  horizontally 
di  rected  . 

The  surface  patch  technique  [  17  ]  can  be  employed  to  solve  this 
boundary  value  problem.  An  electric  field  integral  equation  based  on  the 
boundary  condition 

rix  (E  +  E  )  =  0  on  the  microstrip  patch  (1) 

can  be  written.  In  Equation  (1),  E*-  is  the  field  due  to  the  excitation 

probe,  h  is  a  unit  normal  to  the  patch,  and  Es  is  the  electric  field  due  to 
the  equivalent  surface  current  Js . 

It  is  necessary  to  point  out  that  the  field  at  the  source  point 
contains  a  singular  integrand,  as  can  be  seen  in  the  field  equations  in 
Reference  15  used  to  compute  the  tangential  field.  It  appears  that  the 
singularity  is  of  the  1/ to  type.  No  attempt  was  made  to  deal  with  this 
problem  directly.  Instead,  the  testing  is  performed  slightly  away  from  the 
source  point  which  is  the  technique  used  by  Richmond  [18  )  and  Harrington 
fl  9j  in  wire  antenna  analyses. 

Another  difficulty  is  due  to  the  fact  that  the  source  and  field  points 
have  the  same  z  coordinates.  This  can  be  dealt  with  by  choosing  a  small 
yet  numerically  acceptable  difference  between  the  z  coordinates  of  the 
source  and  the  field. 
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This  type  of  point-matching  method  was  applied  to  microstrip  problems 
by  Farrar  and  Adams  l  20  1  to  calculate  the  capacitances  of  a  patch  and  the 
gap  between  two  patches.  Their  work  was  for  a  much  simpler  static  model 
with  little  direct  relevance  to  the  antenna  problem,  which  requires  a  full 
solution  of  the  time-harmonic  fields.  However,  examination  of  their  work 
in  comparison  with  the  present  approach  indicated  that  development  of  an 
algorithm  in  the  present  approach  would  be  impractical  and  difficult.  In 
addition,  their  work,  as  well  as  others'  l  21-23),  led  to  an  awareness  that 
the  peripheral  magnetic  current  plays  a  predominant  role  in  the  physics 
involved  in  microstrip  structures.  While  earlier  work  addresses  only 
circuit  (interior)  problems,  the  method  of  Okoshi  1  21  ]  has  direct 
implication  to  antenna  problems.  The  equivalent  electric  current  sheet 
method  was  thus  abandoned  in  its  early  stage  of  development, 

C.  Solution  in  Terms  of  Peripheral  Magnetic  Current 

The  majority  of  analyses  on  microstrip  structures  take  advantage  of 
the  fact  that  the  substrates  are  electrically  thin.  As  a  result,  the 
fields  in  the  region  between  the  patch  and  the  ground  plane  contain 
primarily  vertical  E  (perpendicular  to  the  ground  plane)  and  horizontal  H 
fields.  As  shown  in  Figure  3,  the  fields  in  this  region  are  also  quite 
independent  of  the  z  coordinate.  At  the  edge  of  the  microstrip,  current  on 
the  patch  normal  to  the  edge  must  vanish  as  required  by  the  continuity  of 
current.  As  a  result,  the  tangential  component  of  H  field  along  the  edge 
is  negligible.  Thus,  the  peripheral  slot  can  be  represented  by  an 
equivalent  magnetic  current  n  x  Es ,  as  shown  in  Figure  3,  where  n  is  a  unit 
vector  normal  to  the  peripheral  surfaces  (n  is  always  perpendicular  to  z), 
and  Es  is  the  total  electric  field  on  s  given  by 

Is  =  Ez  z  .  (2) 


In  this  approach,  Eg, or  the  equivalent  magnetic  current  n  x  Es ,  is  the 
crux  of  the  problem.  This  model  has  been  extensively  reported  in  the 
literature  in  conjunction  with  other  numerical  techniques  for  microstrip 
circuit  problems.  For  microstrip  antennas,  however,  it  was  primarily  used 
in  the  approximate  computation  of  the  far-zone  fields.  It  is,  therefore,  a 


natural  step  to  expand  this  technique,  which  has  been  very  successful  in 
dealing  with  microstrip  circuits,  to  microstrip  antennas. 

The  approach  of  Okoshi  and  Miyoshi  [21]  appeared  to  be  particularly 
suitable  to  antenna  problems  because  it  deals  directly  with  the  exterior 
problem  with  an  integral  equation.  Consequently,  Gupta  visualized  the 
advantage  of  this  approach  and  applied  it  to  microstrip  antennas  I  24-27  I. 

An  additional  advantage  of  the  peripheral  magnetic  current  method  is 
its  adaptability  to  the  microstrip  antennas  with  radial  conducting  strips 
for  dual  frequency  operation  invented  by  Mcllvenna  and  Kern'ois  [10,11  I. 
Effects  of  the  radial  conducting  strips  can  be  treated  as  a  circuit  problem 
in  which  the  impedance  of  the  strip  is  computed  with  closed-form  analytical 
formulas . 

A  new  approach  was  taken  in  this  research.  The  exterior  field  of  a 
microstrip  antenna  is  represented  by  a  magnetic  current  over  he  peripheral 
opening  and  a  conducting  surface  enclosing  the  antenna  region  as  shown  in 
Figure  4,  [2,13  1  .  This  model  follows  directly  from  the  equivalence 

principle  and  the  uniqueness  principle  (db,  p.  108  [. 

As  shown  in  Figure  5,  we  can  now  expand  Ms  into  a  peripheral 
distribution  along  s,  the  periphery,  as 

N 

Ms  =  2  ^snP(Sn),  (3) 

n=l 

where  P(Sn)  is  a  pulse  function  such  that, 

P(Sn)  =  1  if  Sn>S>Sn+1,  (4) 

=  0  elsewhere  . 

Msnis  a  horizontal  magnetic  current  element  that  produces  fields  as  if 
the  microstrip  antenna  region  did  not  exist.  That  is,  Msn  produces  fields 
in  the  substrate  over  the  ground  plane  --  a  problem  treated  extensively  in 
the  preceding  section. 

The  presence  of  the  microstrip  antenna  region,  represented  by  a 
conducting  pill  box,  in  Figure  4  requires  that  the  amplitude  and  phase  of 
all  the  Msn's  be  adjusted  so  that  the  tangential  E  fields  vanish  on  the  top 
and  the  side  (periphery)  of  the  pill  box.  Note  that  fields  due  to  each  M,.  _ 


Segmentation  of  the  periphery  of  a  microstrip  antenn 
of  arbitrary  geometry. 


satisfy,  by  definition,  Maxwell's  equations  and  the  boundary  conditions 
outside  (but  not  on)  the  pill  box. 

Note  that  in  this  case  it  is  not  necessary  to  enforce  the  boundary 
condition  for  tangential  H  even  though  in  general  a  boundary  condition  fni 
tangential  H  is  also  required.  This  is  due  to  the  fact  that  in  the  present 
formulation,  tangential  H  on  the  dielectric  interface  due  to  each  Msn  is  by 
definition  already  continuous  at  the  interface.  Thus  an  integral  equation 
of  the  form 


ill  Ms  +  Mi  1=0  (5) 

on  the  microstrip  patch  can  be  written.  is  the  end-fed  excitation 

magnetic  current  of  the  microstrip  antenna. 

Physically,  this  formulation  is  similar  to  those  of  Okoshi  and  Gupta 
and  their  colleagues  121,24-27].  The  integral  equation  is  expressed  in 
terms  of  the  peripheral  magnetic  current  and  excitation  current  alone. 
This  is  a  definite  advantage  in  treating  arbitrarily  shaped  patches  and  in 
particular,  the  antenna  with  radial  conducting  strips. 

The  present  method  can  also  be  viewed  as  a  modal  or  eigenfunction 
solution.  All  the  possible  modes  of  fields  for  the  exterior  and  interior 
regions  are  those  due  to  horizontal  magnetic  current  elements  at  the 
periphery.  The  combined  effect  of  the  microstrip  patch  is  that  tangential 
K  field  vanishes  on  this  microstrip  patch  of  infinitesimal  thickness.  This 
approach  is  similar  to  the  modal  expansion  techniques  used  extensively  in 
waveguide  discontinuities  and  array  antennas. 

A  computer  program  based  on  Equations  (A)  and  (5)  was  written  to  solve 
the  microstrip  antenna  problem  by  the  method  of  moments  l  19  ].  However,  a 
lack  of  anticipated  numerical  convergence  was  observed.  During  the 
exhaustive  debugging  stage,  it  became  clear  that  new  concepts  and 
techniques  were  involved  in  this  computer  program  which  needed  more 
debugging  and  investigation.  This  promising  method  in  dealing  with 
none  anon i c a  1  microstrip  antennas  merits  further  extensive  numerical 


SECTION  III 


ANALYSIS  OF  A  MICROSTRIP  ARRAY 


A.  Introduction 

The  analysis  of  a  microstrip  array  is  a  difficult  problem.  As  a 
result,  very  little  has  been  reported  in  this  area  even  though  other 

arrays  such  as  waveguide  and  dipole  arrays  have  been  extensively  studied. 

Thorough  reviews  on  microstrip  arrays  can  be  found  in  references  2  and  29. 
The  difficulties  in  microstrip  array  analyses  are  primarily  in  the  handling 
of  the  interaction  in  the  feed  region.  Recently,  Liu,  et  al  [  3 0 1  reported 

an  analytical  method  which  deals  with  the  feed  by  approximating  the 

vertical  probe  feeds  with  a  continuous  and  uniform  current  sheet. 
Interestingly,  they  also  broke  up  the  analysis  into  two  simpler  ones  by 

employing  the  equivalence  principle  —  a  technique  similar  in  many  aspects 
to  the  peripheral  magnetic  current  method  described  in  the  preceding 

section.  Since  future  development  of  phased  arrays  will  probably  emphasize 
monolithic  and  multilayer  structures  1 2 9 1  * emphas is  in  this  research  was  on 
the  PIN  diode  switched  structure  shown  in  Figures  6  and  7.  The  tinv  gaps, 
each  about  10  mils  wide,  are  connected  with  PIN  diodes  which  are  switched 
on  and  off  at  a  low  rate.  The  long  thin  slots  in  Figure  6  insulate  the 
diode  bias  currents  which  flow  through  the  conducting  (black)  ..urfaces. 
The  slots  (in  white)  are  switched  on  and  off  with  these  diodes. 
Alternatively,  conducting  strips,  instead  of  slots,  can  be  designed  to  form 
a  switching  surface,  as  shown  in  Figure  7.  It  is  expected  that  the 

techniques  developed  in  treating  these  types  of  microstrip  arrays  will  be 
useful  for  other  similar  microstrip  arrays.  Thus,  a  more  general  type  of 
element  slot  configuration  as  shown  in  Figure  8  was  chosen  for  the  computer 
program  input.  The  switching  of  the  phase  front  for  the  reflected  or 
transmitted  wave,  and  thus  the  direction  of  the  antenna  beam,  is 

accomplished  by  the  individually  switched  PIN  diodes. 

There  are  two  major  analytic  difficulties  involved  in  the  computation 
of  scattering  of  this  configuration;  the  complex  geometry  of  the  aperture 
and  the  effect  of  the  diodes.  Scattering  of  screen  structures  for 
rectangular  and  circular  apertures  (or  plates)  were  treated  by  Chen  [  31  ,  3 2 1 
and  Lee  [ 3  3 i  .  Analyses  involving  loading  and  other  geometries  were 
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conducted  by  Munk,  et  al  [  35  ]  for  thin  slots  and  dipoles.  Recent 

research  in  the  Jerusalem-cross  geometry  has  also  been  performed  for  thin 
crosses  for  which  simple  current  distributions  can  be  reasonably  postdated 
l ib-  id  ]  . 

In  this  research,  the  effects  of  loading  and  slot  geometry  are 
included  —  the  slot  is  not  assumed  to  be  thin.  It  is  difficult  to  find  a 
suitable  set  of  basis  functions  to  represent  the  unknown  aperture  field  to 
obtain  a  numerical  solution.  Natural  waveguide  modes  for  this  slot  could 
be  obtained  in  a  manner  similar  to  the  ridged  waveguide  solutions  by 
Montgomery  [id  ].  However,  the  analysis  would  be  very  complex  ind  this  is 
why  it  has  not  yet  been  attempted.  Another  possible  approach,  which  was 
chosen  in  this  research,  is  to  expand  the  unknown  aperture  fields  in  a  set 
of  orthogonal  exponential  basis  functions.  These  exponential  functions  are 
closely  related  to  the  sinusoidal  functions  representing  rectangular 
waveguide  modes,  but  they  are  not  as  restrictive  on  the  boundaries  so  that 
continuity  at  X  =  A/ 2  -  da  is  possible.  However,  difficulties  are  expected 
in  this  method  in  establishing  criteria  for  relative  convergence  [3‘! ,  -'tO  ]. 


£.  B.  Formulation  of  the  Approach 

We  begin  by  dividing  the  problem  into  the  TE  (transverse  electric)  and 
TM  (transverse  magnetic)  cases.  The  incident  wave  for  each  ca-e  can  be 
represented  as 

E*  =  $  (- j  kQs  in0  )  A^  e-j*So‘!  TE  ,  (p) 

and 


E1  =  Q  (jkosin0)  ru  A*  e  j]So*£ 

— m  J  u  m 


TM  , 


wlie  re 


0Q  k u  0/t0  » 

kn  -  k0  [  sinB  cos.})  x  +  s  i  n0  s  i  n(}>  y  +  cos(3  ?.  ]  > 


K>  ~  oh,) 


r  =  rr-  xx  *■  vy  +  zz  , 

(  A  A  A 

U  :  cosy  cos.})  x  +  cos'd  s  i  n>{)  y  -  s  i  nB  z  ,  and 
'1  ~  -siiv?  x  *•  cost})  y  . 
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It  is  possible  to  represent  fields  in  all  four  regions,  I  through  IV 
in  Figure  d  ,  in  terms  of  vector  potential  A  in  each  region  [ill  .  For  the 
incident  wave. 


A1  =  z  A1  =  z  A1  e_j*i’£  ,  and 
— e  e  z  e 


(q  i  ) 


.1  ,  i  ~ 

A  =  z  A  =  z 
— m  mz 


4  e-jfc’I  . 


(9b) 


We  can  write,  for  each  region,  the  forward  and  backward  propagating  vector 
potentials  as  follows. 
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where  the  summations  over  p  and  q  are  actually  truncated  to  ±P  and  4Q,  and 

^pq  (x  >y  )  =  (AB)_l5  e“j(upox  +  vpqy^  > 

upo  =  2np/A  +  k0  sin!  cos,i  ,  and 

Vpq  =  2  'i  q  / 15  +  k0  sinO  sim|) 

A  pair  of  integral  equations  can  be  obtained  by  enforcing  boundary 
conditions  at  the  three  dielectric  interfaces.  Let  the  electric  field  at 
the  aperture  region  at  z  =  0  be 


(  1 4 .  i ) 

(14b) 

(I  4c) 


Hap  E.ax  x  +  ^ay  V  » 


(14) 
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The  following  two  integral  equations  can  be  obtained  by  matching 
field  at  the  boundary  at  z  =0, 


the  E 


yVa y  '^pq  (x»y)  dxdy  =  -jupo  [  5p  6 

a  per 


o  r  o  (1)  .  i  (2)  r 
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q  pq  e  pq  epq 


+  jvpqapq  [  6°  5°  3(n  a1  +  cj(2)  Ar  1  , 
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E  '{>*  (x,y)  dxdy  =  _jVD0  [  5°  5°  A^  A 

ix  ’pq  3  P°  l  P  q  pq  e  pq 


r  1 

epq  J 
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>  jupqapq  f  <S°  A°  n  ^  A*  +  cf^  \T 

—^7 —  l  p  q  m  p^  mP^  (17) 


for  p ,  q  =  0 ,  *1 ,  *2,  ...  (± ,  ±Q)  , 

where  the  integration  is  over  the  aperture  of  a  unit  cell  at  z  =  0  and  1,0 

P 

and  S°  arc  the  Kronocker  delta  function  and  the  coefficients  <i  ,  etc.  are 

q  pq 

defined  in  the  appendix.  Although  the  derivation  so  far  has  been  along  the 
same  path  as  that  of  Lee  [33J,  there  are  nontrivial  differences  in  the 
mathematical  expressions  derived.  Thus,  the  detailed  expressions  are 
included  in  Appendix  B  for  comparison. 

Equations  (16)  and  (17)  contain  too  many  unknowns.  They  can  be 
further  reduced  by  enforcing  the  boundary  condition  for  H  field  at  z  =  0 
and  eliminating  all  the  unknown  A  parameters  according  to  the  transmission 
relationships  in  the  appendix.  Thus,  we  have 
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I,,  explicit  description  of  this  integral  equation  is  included  in  the 
;  pond i x . 
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C.  Moment  Method  Solution 

The  pair  of  integral  Equation  (lH)  can  be  solved  by  the  well  known 
method  of  moments.  Two  major  difficulties  are  involved:  the  construction 

of  suitable  basis  functions  and  the  treatment  of  the  relative  convergence 
problem.  We  choose,  for  the  unknown  aperture  fields,  the  following 
expans  ion: 

M  N 

Eo„  =  7  7  A  ■[’  ( x ,y )  in  region  C  , 

ax  4-4  4m*  mn  ran  ”  b 

m=-M  n=-N 

=  0 
M  N 

E  =  V"'  ^  '  B  <t>  (x,y) 

ay  ZmJ  4mmi  mn  mn 

m=-M  n=-N 

=  E8  =  E  e  j  ( sx+w)  cos  {[|x|  _  (A/2  -  da)]  71/ ( 2g)>  in  D. 
ay  8  (JO) 

where  C  and  D  are  regions  of  the  aperture  in  the  unit  cell  shown  in  Figure 
9  .  g  is  the  effective  gap  width  and 

g  =  gon  when  the  diode  is  forward-biased  »  and 

=  g0ff  when  the  diode  is  reverse-biased  •  (Ji) 

The  basis  function  mn  constitute  an  orthonorraal  set  in  region  C  if  we  choose 

■  mn  <x.y>  =  l  (A  -  2  da)  (B  -  2  db)j-'s  e-j(U'pox  +  V^qy)  ,  (  >2) 

where 

»p0  ~  2pT/(A-2  da)  +  k0  sin90  cos4>0  ,  and  (2  1) 

Vpq  =  -  2  dh)  +  k0  sin«0  siivt>0  . 

The  abo v.  •■xpansion  can  also  be  set  up  in  the  following  manner  to  show  the 

orthogonality  of  t  lie  basis  functions.  We  write 


in  region  D  ,  <  1  •)  ) 

in  C  , 
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where  , 


>!>  =  i>  Pr  ,  and 

mn  mn 


(2b  ) 


y  =  p„  +  6  n  6 N+  *  e®  P 


mn  mn 


m  n  ay  d 


(27) 


Pc  and  Pj  are  pulse  functions  defined  as. 


Pc  =  1 


in  c  , 


(28) 


=  0 


elsewhere . 


X  y  .  .  m 

The  orthogonality  among  T  ^  and  is  readily  obvious  since  Pc  Pj  = 

0.  This  is  the  subsectional  expansion  discussed  in  Reference  19  . 

In  addition  to  the  integral  Equation  (18),  we  need  two  more  equations 
to  insure  the  continuity  of  tangential  E  fields  at  x  =  *  (A/2-da)  and 

|y|<dg/2.  Thus  we  add  the  following  equations 

M  N 

X!E  1  -  lA/2_'‘'’  >'  *  Eay(‘[A/!-d»].  0)  <J'0 

iii~ -M  n  = -N 

The  use  of  exponential  functions  to  represent  fields  in  Region  c  is 
advantageous  in  tile  present  problem  because  it  allows  non -vani shi ng  Ey  at  x 


-  £(A/2-da)  needed  for  the  continuity  of  fields  between  regions  C  and  D. 
However,  the  choice  of  the  exponential  function  is  not  as  efficient  as 
waveguide  modes  in  meeting  the  edge  conditions  at  other  conducting  edges 
away  from  the  gap. 

Solution  of  Equations  (IB)  and  (  29)  was  carried  out  by  using  the  same 
relative  convergence  criteria  established  by  Lee  [  11  ]  for  rectangular 
apertures.  Thus,  P  and  Q  are  first  chosen  based  on  an  initial  guess  of  the 
physical  nature  of  the  problem  and  the  memory  size  of  the  computer.  M  and 
N  are  then  chosen  as 


,  and 


where  [  a  ]  denotes  the  largest  integer  in  a.  As  was  pointed  out  in 

References  19  and  t(),  this  choice  of  M  and  N  for  the  particular  P  and  Q 
leads  to  the  most  accurate  results;  larger  M  and  N  often  yield  poor 
results. 

The  effect  of  the  gap  is  a  reduction  in  the  effective  da  in  Equation 
(  10)  .  For  arrays  with  small  da,  as  often  encountered  in  practice,  the  gap 
has  little  effect  on  Mc ;  Mc  can  be  chosen  according  to  Equation  ('<•)  or 
with  one  additional  mode. 

D  .  Numerical  Com put  at i o n 

A  computer  program  was  written  to  analyze  the  type  of  microstrip  array 
shown  in  Figure  H  with  emphasis  on  the  strip  type  shown  in  Figure  , 

which  is  the  special  case  with  d|,  =  0.  Testing  and  debugging  of  the 

program  began  with  cases  in  which  the  metal  screen  was  absent,  that  is, 

■  l,i  dg  -  0.  The  problem  is  the  wl  1-i.nown  reflection  and  transmission 

through  a  dielectric  layer.  Table  1  shows  a  comparison  between  this 
coriiput  at  i  on  and  the  results  of  Lee  (estimated  from  Figure  2  of  Reference 
1  >)•  The  agreement  can  he  considered  extremely  good  since  Lee's  data  were 
extracted  from  a  figure  by  estimation.  Note  that  there  is  a  IS"' 

difference  between  the  phases  of  TK  and  TM  reflection  coefficient.  This  i  ■; 


TABLE  I 


COMPARISON  FOR  THE  REFLECTION  COEFFICIENT  COMPUTATION  FOR  THE 
CASE  OF  da  =  db  =  dg  =  0  (SCREEN  ALSO)  AND  A  =  0.2,  B  =  0, 

=  1 2  =  0.1  (ALL  IN  WAVELENGTHS)  erj  =  ^r2  =  2-5 


ANGLE 


Tt 


TM 


(deg) 


0 

B 

Amp. 

phase 

Amp. 

m 

M 

(deg) 

III 

1° 

0° 

0.395 

-130 

0.398 

50 

(0.395  ) 

(  -130  ) 

(  0.395  ) 

(  50  ) 

0.54 

-137 

0.267 

42 

o 

o 

0° 

(0.54  ) 

(  -1)7  ) 

(  0.267  ) 

(  42  ) 

due  to  the  difference  in  sign  in  the  incidence  wave  in  Equations  (t>  )  and 
( /'  ).  To  be  a  true  reflection  coefficient,  180°  has  to  be  added  to  the  TM 
case . 

Next  we  check  the  case  in  which  dj,  =  0  so  that  the  screen  becomes  an 
array  of  parallel  strips.  The  dielectric  layers  are  chosen  to  be  absent  by 
letting  er i  =  t  r  2  =  1.  The  computed  results  are  again  compared  with  those 
of  Lee  in  Table  II  with  excellent  agreement.  However,  we  noticed  that  the 
choice  of  M  and  N  is  essential  to  the  accuracy  of  the  results.  For 
example,  if  we  increase  the  number  of  M  modes  by  choosing  M  =  Mc  +  1,  (Mc  = 
1),  the  computed  TM  reflection  coefficient  becomes  0.122  at  v  =  1°  and 

0.0939  at  0  =  40°,  even  though  other  computed  results  change  very  little. 

Thus,  a  14%  change  in  computed  results  is  highly  possible  if  the  choice  of 
modes  is  only  slightly  different.  This  is  the  well-known  phenomenon  of 
relative  convergence.  In  the  following,  it  will  be  shown  that  this  remains 
to  be  a  difficulty  in  the  use  of  this  computer  program. 

Next,  we  let  sr  =  2.5  in  the  two  dielectric  layers  in  the  case  of 

Table  II.  The  results  are  compared  again  with  Lee’s  in  Table  III.  Large 
discrepancies  now  appear,  especially  for  the  TE  case.  To  examine  whether 
this  is  due  to  a  lack  of  relative  convergence,  we  perform  a  convergence 
test  and  present  the  results  in  Figures  10  and  II  for  the  amplitude  and 
phase  of  the  reflection  coefficient.  As  can  be  seen,  the  convergence  is 
fairly  good  near  Mc . 

Recalling  that  the  excellent  agreement  in  Table  II  had  not  been 

reached  until  two  printing  errors  in  the  input  data  described  in  the  text 

of  the  paper  were  detected,  it  was  conjectured  that  errors  might  exist 

also  in  tire  input  data  for  the  case  in  Table  III.  The  discrepancy  was  thus 
believed  to  be  due  to  differences  in  input  data,  not  program  bugs.  To 

further  illustrate  this  point,  let  us  examine  the  case  of  a  full  screen  in 
which  A=B=0.4,  2da  =  2dp=0.2  A,  t(=t2=0.1,  dg=0  (all  in  wavelengths),  and 
It  =  _|r  =  2.5.  The  total  transmission  coefficient  (including  the 

depolarized  wave  transmitted)  in  dB  is  computed  in  comparison  with  Lee’s 
data  in  Figure  l .  The  incident  wave  is  transverse  magnetic  with  its 

incidence  angles  indicated.  The  agreement  is  quite  good  except  at  near 
-58*’,  where  the  present  calculation  has  a  siiarper  drop  than  tiiat  of  Lee. 
An  even  more  challenging  case  is  shown  in  Figure  1  ' ,  in  which  the 
computations  of  the  amplitude  and  phase  of  the  TM  transmission  coefficient 

2f> 


TABLE  II 


COMPARISON  FOR  THE  REFLECTION  COEFFICIENT  COMPUTATION 
FOR  THE  CASE  OF  A  =  B  =  0.2,  dfc  =  dg  =  0,  da  =  0.05, 
tt  =  1 2  =0.1  (ALL  IN  WAVELENGTHS),  Lr i  =  ' r2  =  1 


ANGLE 

(deg) 

TE 

TM 

8 

■ 

Amp. 

wm 

Amp. 

phase 

H 

Issfl 

(deg) 

0.978 

171 

0.142 

82 

1° 

0° 

(  0.9  78  ) 

(  169  ) 

( 0.142  ) 

(  83  ) 

0.992 

173 

0.110 

84 

40° 

0°* 

( 0.990  ) 

(  173  ) 

(o.uz  ) 

(  84  ) 

(  ):  LEE  (estimated) 
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TABLE  III 


COMPARISON  FOR  THE  REFLECTION  COEFFICIENT  COMPUTATION 
FOR  THE  CASE  OF  A  *  B  =  0. 2,  db  =  dg  =  0,  da  =  0.5, 
t 1  =  1 2  =  0.1  (ALL  IN  WAVELENGTHS) ,  fr  [  =  e  r2  =  2- 5 


ANGLE 

(deg) 

TE 

TM 

0 

■ 

Amp. 

phase 

Amp. 

phase 

H 
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art?  compared.  The  sharp  fall  and  rise  of  the  transmission  coefficient  near 
0=60.6°  and  ,£=4 5°  occurs  in  both  computations.  Even  though  the  agreement 
between  0=63°  and  70°  is  not  quite  satisfactory,  the  excellent  agreement 
exhibited  in  this  important  blind-spot  phenomenon  suggests  that  the 
discrepancy  is  probably  not  due  to  program  bugs. 

Although  some  of  the  discrepancies  have  not  been  completely  resolved, 
the  numerical  check  is  considered  satisfactory.  Numerical  testing  for  the 
case  of  dg  =  0  was  handicapped  by  a  lack  of  existing  data  to  compare  with. 
Without  other  data  to  guide  the  initial  analysis,  it  is  difficult  to 
perform  conclusive  numerical  checks.  Indirect  numerical  checks  were 
conducted,  but  they  were  inefficient,  costly,  and  inconclusive  and  are  not 
discussed  in  this  report. 
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Techniques  }  r  t  t  -  in.ilv;  is  of  noncanonical  microstrip  antennas  and 
array;  were  i  nv.-s  1 1  rated  w,t".  the  ultimate  goal  of  antenna  performance 
improvements.  There  were  three  major  thrusts  in  this  research,  each  of 
which  culminated  in  the  development  of  a  particular  computer  program. 

The  general  problem  of  the  radiation  of  an  arbitrary  Hertzian  dipole, 
whether  horizontal  or  vertical,  electric  or  magnetic,  in  general  stratified 
media  was  numerically  solved.  The  media  can  have  multilayers  of  various 
thicknesses  with,  complex  dielectric  constants.  The  shortcomings  of  this 
approach  are  that  it  is  not  directly  applicable  to  far  fields  and  its 
convergence  depends  on  the  vertical  (z)  spacing  between  the  dipole  and  the 
field  points. 

The  Hertzian  dipole  algorithm  was  developed  as  a  basis  for  the 
analysis  of  a  general  noncanonical  microstrip  antenna.  A  method  was 
developed  to  treat  the  antenna  solely  in  terms  of  a  horizontal  magnetic 
current  along  the  periphery  of  the  patch.  This  moment  method  has  several 
advantages,  especially  in  dealing  with  microstrip  antennas  with  radial 
conducting  strips.  The  computer  program  development  encountered 

difficulties  in  numerical  convergence.  ,  Because  several  new  concepts  and 
techniques  are  involved  in  this  approach,  further  extensive  examinations 
are  needed. 

Microstrip  array  analyses  were  focused  on  a  monolithic  di ode-swi t ched 
screen  sandwiched  in  two  dielectric  layers.  A  moment  method  approach  with 
exponential  basis  functions  was  taken.  The  diode  gap  region  is  assumed  to 
have  a  single  mode  sinusoidal  field  distribution.  The  computer  program  was 
tested  against  data  in  the  literature.  Although  some  excellent  agreements 
were  observed,  discrepancies  also  appeared  which  need  further 
i nves t i gat i on . 

This  research  program  is  an  ambitious  endeavor  as  it  deals  with 
difficult  electromagnetic  and  numerical  analysis  problems.  The  present 
accomplishments  are  expected  to  contribute  to  integrated  circuits  and 
geophysical  science,  in  addition  to  microstrip  antennas  and  arrays.  For 
all  three  major  tasks,  continued  research  is  recommended  for  full 
development  and  validation  of  these  new  techniques. 
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Relationships  among  the  vector  potentials  are  : 
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Propagation  and  dielectric  parameters  are,: 
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